Quasi-normal modes for black hole solutions unknown in an analytical form 
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We review the papers 0, Q, 0J- We discuss possibilities of studying the quasi- normal modes of 
black holes that are not known in an analytical form. Such black holes appear as solutions in various 
theoretical models and real astrophysical approximations when one takes into account the black hole 
neighborhood. 
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I. INTRODUCTION 

There are three stages of the perturbation evolution 
near a black hole. The initial black hole outburst from 
the perturbation source takes place at early time of 
the evolution. Then, at late time, the damping (quasi- 
normal) oscillations appear, finally suppressed by power- 
law or exponential tails at asymptotically late time. The 
quasi-normal modes (QNMs) can be characterized by a 
set of complex frequencies: real part is the actual oscilla- 
tion frequency while imaginary part describes the damp- 
ing rate of the oscillation j4j. 

The QNMs do not depend on initial perturbations be- 
ing, thereby, an important characteristic of a black hole. 
They are expected to be observed in the nearest future 
with the help of a new generation of gravitational anten- 
nas. Therefore, it is important to know how they change 
when the black hole solution has some corrections due to 
some matter and fields or due to the possible existence 
of some physics beyond the classical Einstein gravitation. 
Mostly, such solutions could be obtained only numeri- 
cally by resolving some system of differential equations. 

According to AdS/CFT (anti-de Sitter/Conformal 
Field Theory) correspondence, QNMs of a large black 
hole in AdS space-time should coincide with the poles of 
the retarded Green function of the perturbation in Af = 4 
SU(N) super- Yang-Mills theory for large N That is 
why asymptotically AdS black hole solutions are being 
extensively studied as well. 

It turns out that the dominating frequencies depend 
mostly on the black hole solution behavior in some re- 
gion near the event horizon. Thus one can find them in 
frequency domain even if the large distance behavior of 
the solution is not known at all. 

We consider two qualitatively different examples: hairy 
black hole in anti-de Sitter background and asymptoti- 
cally flat Einstein- Aether black hole solution. We will 
show that in both cases the behavior of the solution at 
large distances is not important because of different rea- 
sons: 

• For asymptotically anti-de Sitter background we re- 



quire Dirichlet boundary conditions at spatial infin- 
ity. The most significant part of the metric pertur- 
bations stays, thereby, near the black hole. That 
is why the solution behavior at this region causes 
dominant influence on the QN spectrum. 

• For the asymptotically flat case the QNMs search- 
ing can be reduced to the scattering problem. Be- 
ing resonant oscillations that survive at late time, 
the QNMs in this approach are the poles of the re- 
flection coefficient. Therefore, the QN frequencies 
are determined mainly by the form of the effective 
potential near its peak. 



II. QN SPECTRUM OF A BLACK HOLE IN 
ADS SPACE-TIME 

Since our approach for the QN modes calculation of 
asymptotically AdS black hole solutions is based on the 
Horowitz-Hubeny method Q let us start from a brief 
review of it. 

Suppose we have an analytical solution for a spher- 
ically symmetric black hole metric. Then any spheri- 
cally symmetric perturbation after separation of variables 
$(t, r) = exp(— iu))(j)(r) satisfies a second order differen- 
tial equation which has the form 



<j)"(r) + A(r)(f>'{r) + B(r)4>{r) = 0. 



(1) 



One should note, that in many cases the equation for 
the radial part can be reduced to the form {1} for the 
non-spherically symmetric perturbations as well. 

By definition, QNMs are eigenvalues w of |l} with 
boundary conditions (BCs) which require ingoing wave 
at the event horizon and zero at spatial infinity [5[ . Such 
BCs allow us to exclude a regular singularity of the equa- 
tion |T|) at the event horizon. After introducing the new 
variable z, that is zero at event horizon and one at spatial 
infinity, one can obtain the equation 



s(z)z 2 y"(z) + t(z)zy'(z) + u(z)y(z) 



0. 



(2) 
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where y(z) = <t>(z)z tuJK and n > is chosen to make y(z) 
regular at even horizon (z = 0). 

The Horowitz-Hubeny method works if s(z), t(z) and 
u[z) are polynoms and all the other singularities of the 
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equation lay outside the unit circle [1C 
can expand y{z) as 



In this case one 



l{z) = ^2y n z n , y = l. 



(3) 



n=0 



The coefficients of the series satisfy the recurrence rela- 
tion 



n-l 



yk(k(k - l)s n - k + ktn-k + Un-k) ^ 



fc=0 



n.(n — 1)sq + nto 



which is finite since s(z), t(z) and u(z) are polynoms: 



s(z) 



E< 

71=0 



t(z) 



E 

n=0 



u(z) = 



Thus one can find w by minimizing |y(l)| restricted by 
some large number of expansion terms in (|H). 

It is clear that if s(z), t{z) and u(z) are series and 
s(l), t(l) and w(l) converge quickly enough, then y n is 
still possible to calculate because it depends insignifi- 
cantly on higher terms of the series. In order to find 
series expansion for s(z), t(z) and u(z), one can use the 
equations which define the black hole solution. We can 
always do this because s(z), t{z) and u(z) can be explic- 
itly expressed in terms of the metric coefficients and their 
derivatives [lj. 

It is important to note that one can control the preci- 
sion of the eigenvalues uj by requiring the convergence of 
the found result with respect to increasing of the number 
of expansion terms of all the series. 

This technique was used in [l[ to study QNMs of the 
scalar hairy black hole in the AdS background @ . 



III. QN SPECTRUM OF A BLACK HOLE IN 
THE FLAT BACKGROUND 



For asymptotically flat black hole solutions, the pertur- 
bation equation is usually reduced to a wave- like equation 



-) = 0. 



(5) 



The general solution of the equation {5| at infinity is 

* = Ai n tpi n + Aoutipout, r* — > 00. 



(6) 



The quasi-normal modes in this approach, by the defini- 
tion, are the poles of the reflection coefficient A out /Ai n . 
Our starting point is suggested by the WKB method [7[ 
where the asymptotic solutions of the wave equation near 
the event horizon and near spatial infinity are matched 
with Taylor expansion near the peak of the potential, i.e. 
between the two turning points V(r) — lo 2 = 0. 

We conclude therefore that the low laying quasi-normal 
modes are determined mainly by the behavior of the ef- 
fective potential near its peak. We shall check that our 



FIG. 1: Potential for electromagnetic perturbations near the 
Schwarzschild black hole (r^ = 1, £ = 2) and the same po- 
tential interpolated numerically near its maximum. Even de- 
spite the behavior of the two potentials are very different in 
the full region of r, except for a small region near black hole, 
low-laying quasi-normal modes for both potentials are very 
close. 




statement is true by considering the well-known potential 
for the Schwarzschild black hole V(x) and also two other 
potentials which lay closely to the Schwarzschild poten- 
tial near its maximum, but has very different behavior 
far from a black hole. These two potentials are chosen in 
the following way. We make a plot of the function V(x) 
of an analytic Schwarzschild potential, then we find some 
number of points for this potential V r near its maximum 
which serve us a basis for our first potential Vi n t which 
is an interpolation of these points near the maximum by 
cubic splines (see Fig. Illlj) . Second potential Vfu is a 
fit of the above plot near the maximum by a ratio of 
polynomial functions. 

The 6-th order WKB formula reads 



iQo 



i=6 

E A > = 

i=2 



f 

2' 



(7) 



where the correction terms of the i-th WKB order A,- 



can be found in [8|, Q = V 



and Qq means the i-th 



derivative of Q at its maximum. For the WKB formula 
of sixth order one needs to calculate derivatives of the 
potential up to 12th order. 

Since the WKB formula contains the value of the ef- 
fective potential its derivatives at the potential peak we 
find that the results obtained with the help of all three 
potentials lay very close if the interpolation and the fit 
were made with the appropriate precision. Despite the 
higher derivatives of our interpolation potential are not 
defined, we are able to evaluate them step by step by in- 
terpolating the first and all the consequent derivatives of 
the potential in the same way. Unfortunately, the inter- 
polation potential is very sensitive to numerical errors. 
Therefore, to calculate the QNMs with the appropriate 
precision, one must find the values of the potential with 
very high accuracy. In fact, for the practical purposes 
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one can use fitting of the potential which does not accu- 
mulate the numerical error. 

To test the accuracy of this approach one can use the 
convergence of the WKB formula 0j. Since the higher 
WKB order depends on the higher derivatives of the ef- 
fective potential, that are more sensitive to the interpo- 
lation or fitting accuracy, the higher order WKB formula 
should give some random values, if the accuracy is not 
enough. 

The possibility of making use of the numerical interpo- 
lation or fitting of the potential near its maximum gives 
us possibility of finding QN spectrum of fields near the 
solutions which are not known in analytical form. This 
technique was used in 0, 0] to find QNMs of the Einstein- 
Aether black hole [§] . The obtained results coincide with 
those calculated in time domain Q . 

Using this technique we are able to make the next step 
in studying of the real astrophysical black holes. Such 
black holes are usually situated in a system of other as- 
trophysical objects (accretion disc, galaxies, dark matter 
etc.). Any such matter causes influence on the metric 
and, therefore, changes QN spectrum of the black hole. 
Since the presented technique does not require neither 
analytical solution nor its behavior at large distance, it 
can be used to estimate the influence of the black hole 
surroundings. 



IV. SUMMARY 

To summarize the above information let us briefly re- 
mind the main points: 

• The low-laying QNMs depend mostly on the black 
hole solution behavior in some region near the event 
horizon. 

• The developed technique allows us to find them if 
the solution is not known analytically. 

• Within the considered approach we are able to es- 
timate the precision of the obtained results. 

• The approach can be used to find QN spectrum for 
a wide class of black hole solutions. 
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